Doping evolution of the absolute value of the London penetration depth and 
superfluid density in single crystals of Ba(Fei_a.Coa,)2As2 
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The zero temperature value of the in-plane London penetration depth, XabiO), has been measured 
in single crystals of Ba(Fei_a;Coa;)2As2 as a function of the Co concentration, x, across both the 
underdoped and overdoped superconducting regions of the phase diagram. For x > 0.047, \ab(Q) 
has been found to have values between 120 ± 50 nm and 300 ± 50 nm. A pronounced increase 
in Aai,(0), to a value as high as 950 ± 50 nm, has been observed for x < 0.047, corresponding to 
the region of the phase diagram where the itinerant antiferromagnetic and superconducting phases 
coexist and compete. Direct determination of the doping-dependent Aa6(0) has allowed us to track 
the evolution of the temperature-dependent superfluid density, from which we infer the development 
of a pronounced superconducting gap anisotropy at the edges of the superconducting dome. 

PACS numbers: 74.25. Nf,74.20.Rp, 74.20. Mn 



INTRODUCTION 

The zero temperature value of the London penetration 
depth is directly related to the superfluid density in the 
ground state of a system through A(0) cx l/y^njfi) 
In the clean, low scattering limit, ns(0) is equal to the to- 
tal density of conduction electrons, n^. There are cases 
in which other phases, for example itinerant magnetism, 
can compete with superconductivity for the same conduc- 
tion electrons, thus reducing the overall number of carri- 
ers in the superconducting state at T = 0. Given the rich 
doping phase diagram of the newly discovered iron-based 
superconductors in which a long range magnetically or- 
dered state, with itinerant character, coexists with a su- 
perconducting state, questions are raised regarding the 
effects of the cornpetition between these states for the 
same electrons [^-Isl ■ One way to approach this matter is 
to study the doping evolution of Aab(O) across the phase 
diagram of these materials and use it to infer the corre- 
sponding change in the superfluid density, especially in 
the regime of the phase diagram where these two phases 
coexist. Determination of the absolute value of the Lon- 
don penetration depth is also important for the correct 
evaluation of the normalized, temperature-dependent su- 
perfluid density, ps(T) = (A(0)/A(r))^. This quantity 
can be calculated from various models of the supercon- 
ducting gap and provides insight into the pairing mech- 
anism. 

In the present study we focus on Aafc(O), which is the 
ground state screening length associated with supercur- 
rent flowing in the crystallographic ab-plane as a result of 
an external magnetic field applied along the c-axis. For 
X > 0.047, the measured values of Aa6(0) have been found 
between 120 ± 50 nm and 300 ± 50 nm. A pronounced 
increase in Aab(O) to a value as high as 950 ± 50 nm 



for X < 0.047 has been observed. We interpret the in- 
crease in Aab(O) for samples with x < 0.047 to be due to 
the competition between the superconducting and itin- 
erant antiferromagnetic states for the same conduction 
electrons. 

The experimental determination of A(0) is a rather 
challenging task since only finite temperatures can be 
reached. There are techniques that are capable of ob- 
taining an estimate of its value by taking advantage of 
the small variation of A(r) at low temperatures, which 
can be on the order of 1 nm/K, along with precision 
measurements. One such technique is muon spin ro- 
tation (/.tSR) [l3|, which has produced estimates for 
AaijfO) of 320 nm in (Bai_^K^)Fe2As2 (T^ ~ 32K) 
[nllll, 470 nm in (Bao.55Ko.45)Fe2As2 (T, ~ 30 K) [H, 
230 nm in Ba(Feo.6Kii4)2As2 (T^ ~ 38 K) [Ti|, 250 nm 
in La(Oi_xFj;)FeAs 115'] and values ranging from 189 nm 



Mini- 



to 438 nm in the Ba(Fei_2;Coa;)2As2 series 

Another technique, magnetic force microscopy, has re- 
ported Aaf,(0) = 325±50nminBa(Feo.95Coo.o5)2As2 [H. 
In addition, optical reflectivity measurements have been 
used to estimate Aaf,(0) in Ba(Fei_xC0:E)2As2 and re- 
ported values of 277±25 nm for x — 0.06 and 315±30 nm 
for for X — 0.08 [19]. It is important to compare the val- 
ues of A(0) obtained by as many different techniques as 
possible because each experiment requires its own set of 
assumptions and modeling procedures. 

Given the overall disparity between the measured val- 
ues of A(0) from these different experimental techniques, 
it is valuable to perform a systematic study of A(0) as 
a function of doping in the series of which large, high 
quality single crystals having homogeneous doping lev- 
els are available, namely the Ba(Fei_2:Co2:)2As2 series. 
The samples used for this study were members of the 
Ba(Fei_a;CO:r)2As2 series and were obtained from the 
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FIG. 1: Scanning electron microscope images of the Al coated 
samples, (a) Large scale view. The broken side is on top. (b) 
and (c) are zoomed in on the Al film on the edge of the broken 
side, (d) A trench produced by a focused-ion beam (FIB), (e) 
Close-up view of the FIB trench showing the Al film and its 
thickness. 



same source as in Ref. [3i] using the same growth pro- 
cedure. 



described in detail in Ref. '^l']. However, as proposed in 
the same reference, the TDR technique can be extended 
to obtain the absolute value of the penetration depth, 
A(r). The key to obtaining A(0) from TDR measure- 
ments is to coat the entire surface of the superconductor 
under study with a thin film of a conventional supercon- 
ductor having a lower critical temperature and a known 
value of A(0). For this study, the aluminum films that 
were used to coat the Ba(Fei_a;Coa;)2As2 samples had 
Tf^ sa 1.2 K and thicknesses of 100 nm, as shown in 

Fig.m 

While the Al film is superconducting it participates 
with the coated superconductor to screen the magnetic 
field generated by the TDR coil. However, when it be- 
comes normal it does effectively no screening because its 
thickness, i, is much less than the normal state skin depth 
at the TDR operating frequency of 14 MHz, where 8ai ~ 
75 /im for p^' = 10 [ifl-cnv [i^l- By measuring the fre- 
quency shift upon warming from Tmim which is the base 
temperature of the sample, to T > T^' we obtain the 
quantity L = XeffiT^'') - \eff{Tmin), shown in Fig. [21 
This quantity can be used to calculate A(0) along with 
the previously determined power-law relation for iron- 
based superconductors [ll], AA(r) — f5T"-, and by using 
the formula for the effective magnetic penetration depth 
into both the Al film and the coated superconductor for 
T < T^'-, which is given by 



EXPERIMENTAL 

The experimental apparatus used for obtaining all of 
the penetration depth measurements in this work was a 
tunnel diode resonator (TDR) J2d\. The essential com- 
ponents of the TDR are a tank circuit formed by an 
inductor and a capacitor, which has a resonance fre- 
quency /o — 1/2tt^/LC w 14 MHz, and a tunnel diode. 
While the diode is biased appropriately it serves as an 
ac power source for the tank circuit. To perform pen- 
etration depth measurements, the sample is mounted 
on a sapphire stage and inserted into the inductor coil. 
The magnetic field of the coil, which is « 10 mOe, is 
screened by the sample and thus changes the induc- 
tance, L, and therefore also the resonance frequency by 
an amount A/. By utilizing Af{T) = -G4nxiT) = 
G[l - (A(T)/i?)tanh(i?/A(T))], the TDR is capable of 
measuring the variation of the penetration depth in a su- 
perconductor, AA(T) — A(T) — A(0), with a resolution of 
nearly 1 A, where G is a geometry dependent calibration 
factor depending on the coil volume, sample volume, de- 
magnetization and empty coil resonance frequency. This 
calibration factor is measured directly by exctracting the 
sample from the inductor coil at its base temperature. 

The TDR technique, as described above, provides very 
precise measurements of the variation of the penetration 
depth, AA(T), but not the absolute value due to reasons 



''"^^^ ^ '^'^'^^.,(T) + A(T)tanh4^' 

where A(T) is the penetration depth of the coated su- 
perconductor and Ayi/(T) is the penetration depth of 
the Al film. As usual with the TDR technique, the 
variation of the penetration depth with temperature, 
AAe//(r) = Ae//(T) - Xeff{Tmin), IS mcasurcd. This 
method has been successfully demonstrated on several 
cuprate superconductors [ilj and has shown agreement 
with measurements of A(0) in FeiJ-j/Tei_a;Sea:) crystals 
obtained by different techniques [24| . Here we use an 
extended analysis obtained by solving the appropriate 
boundary value problem. 

The aluminum film was deposited onto each sample 
while it was suspended from a rotating stage by a fine 
wire in an argon atmosphere of a magnetron sputtering 
system. The formation of non-uniform regions in the 
film was avoided by bonding the wire to only a portion 
of the narrowest edge of each sample. Each film thickness 
was checked using a scanning electron microscope in two 
ways, both of which are shown in Fig.[T] The first method 
involved breaking a coated sample after all measurements 
had been performed to expose its cross section. After 
this, it was mounted on an SEM sample holder using 
silver paste, shown in Fig.[T]Ja). The images of the broken 
edge are shown for two different zoom levels in Fig. [Ijb) 
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FIG. 2: (Color online) Main frame : Full superconducting 
transition of an optimally doped Ba(Feo.g3Coo.o7)2As2 crystal 
before and after coating. Inset : Zoomed in low-temperature 
region, Tm.n < T < before ( green triangles) and after 

(red circles) the Al coating on the same sample. The overall 
frequency shift through the Al transition, denoted as L, is 
used for the calculation of XabiQ)- 



and (c). The second method used a focused ion beam 
(FIB) to make a trench on the surface of a coated sample, 
with the trench depth being much greater than the Al 
coating thickness, shown in Fig. [IJd). The sample was 
then tilted and imaged by the SEM that is built into the 
FIB system, shown in Fig.[TJe). 



RESULTS AND DISCUSSION 

The values of Ao6(0) that were obtained using the pro- 
cedure described above for the Ba(Fei_^CO:E)2As2 sys- 
tem are shown in the top panel of Fig. [3] for doping lev- 
els, X, across the superconducting region of the phase 
diagram, shown schematically in the bottom panel of 
Fig. 131 The size of the error bars for the Aa6(0) points 
was determined by considering the film thickness to be 
i = 100 ± 10 nm and AAi(O) = 50 ± 10 nm. The discrep- 
ancy in Aab(O) for the two samples having x — 0.038 
is clearly beyond these error bars and may possibly 
arise from cracks or inhomogeneities in the Al film, even 
though great care was taken to eliminate them during the 
coating process. Thus, the error bars represent the un- 
certainty of the known parameters and the scatter in the 
data may arise from uncontrolled effects such as cracks 
or inhomogeneities in the aluminum films. The scatter 
in the Aa6(0) values shown in the upper panel of Fig. [3] 
has an approximately constant value of 0.2 /im for all 
values of x, which probably indicates that the source of 
the scatter is the same for all samples. For comparison. 



Fig. |3] also shows Aab(O) obtained from /xSR measure- 
mentsfred stars) [lil . \n\ . the MFM technique (black 
star) [Ig and optical reflectivity (purple stars) [13|, all 
in the Ba(Fei_a;Co2;)2As2 system, most of which are con- 
sistent with our results within the scatter. It may also 
be important to note that the Aah(O) values from other 
experiments are all on the high side of the scatter that 
exists within the TDR AQb(O) data set. This is because 
any cracks or voids in the Al film will lead to underesti- 
mated values of AQb(O). We note that we did not observe 
an increase in Aah (0) towards the overdoped regime as re- 
ported from ^SR measurements [l7| , although our values 
at the optimal doping do agree well. 

Specifically, an increase in Aab(O) on the underdoped 
side below x « 0.047 has been observed, which is in the 
region where the itinerant antiferromagnetic and super- 
conducting phases coexist, as is shown in the bottom 
panel of Fig. [3] The dependence of Aab(O) on carrier con- 
centration is \ab (0) oc 1/ yjnjfi), where is the super- 
fluid density, which is equal to the normal state carrier 
concentration in the clean case. The proportionality be- 
tween Aa6(0) and ns(0) still holds if scattering is included, 
but ris is reduced due to a residual density of states within 
the gap. Overall, an increase in Aaf,(0) is consistent with 
a decrease in the superfluid carrier concentration. There 
is compelling evidence suggesting that the itinerant anti- 
ferromagnetic spin density wave state in these materials 
acts to gap a portion of the Fermi surface 0-0], which 
would remove mobile charge carriers and this qualitative 
idea is consistent with our experimental observations of 
the doping dependence of Aafc(O). Changes in the Hall 
coefficient for these materials, moving from the pure su- 
perconducting region to the coexistence region, have also 
been interpreted as being due to the interaction between 
these phases 2^, 2^ . It has been shown that the opening 
of a superconducting gap in the antiferromagnetic state 
transfers optical spectral weight from a mid-infrared peak 
to a Drude peak, even when the reconstructed Fermi sur- 
face would be fully gapped [13]. As a result, the coexis- 
tence state has a finite n^, although smaller than in the 
pure superconducting state. 

In order to provide a more quantitative explanation 
for the observed increase in Aaf,(0) as x decreases in the 
underdoped region, we have considered the case of su- 
perconductivity coexisting with itinerant antiferromag- 
netism [oj . For the case of particle hole symmetry (nested 
bands), the zero temperature value of the in-plane pene- 
tration depth in the region where the two phases coexist 
is 



'(0) = A°,(0) 




(2) 



where A°j,(0) is the value for a pure superconducting 
system with no magnetism present, and I^af and Ap 
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are the zero temperature values of the antiferromag- 
netic and superconducting gaps, respectively. Deviations 
from particle hole symmetry lead to a smaller increase in 
Aff +'^^^(0), making the result in Eqn.Han upper esti- 
mate. For more information on the details of the calcula- 
tion and the values of A.af and Aq used, see Ref. 0, H^] ■ 
The three blue dashed lines shown in the top panel 
of Fig. [HI which were produced using Eqn. [21 show the 
expected increase in Xab{0) in the region of coexisting 
phases below x « 0.047 by normalizing to three different 
values of Aab(O) in the pure superconducting state, with 
those being 120 nm, 180 nm and 270 nm. This theory 
does not take into account changes in the pure supercon- 
ducting state, so for x > 0.047 the dashed blue lines are 
horizontal. These theoretical curves were produced us- 
ing parameters that agree with the phase diagram in the 
bottom panel of Fig. [3] 0, [1^, which includes a shift of 
the coexistence region to lower values of x by an amount 
of 0.012, and given the simplifications of the model, the 
agreement with the experimental observations is quite 
reasonable. 

While the exact functional form of A(a;) is unknown, 
the solid gray line in Fig. [3] serves as a useful guide to the 
eye (of the form A+B/x"^), which does indeed illustrate 
a dramatic increase of Aab(O) in the coexistence region 
and also a relatively slight change in the pure supercon- 
ducting phase. It should be noted that a dramatic in- 
crease in Aafc(O) below x ~ 0.047 cannot be explained by 
impurity scattering, which would only lead to relatively 
small corrections in the magnitude of Aab(O) (but, indeed, 
significantly affects the temperature dependence of A(T) 

El). 

Values of Aab(O) obtained here can be used to calcu- 
late the actual penetration depth, Xab{T) = AAab(T) -t- 
Aa6(0), where AXab{T) has been measured for each 
Ba(Fei_xCox)2As2 crystal used in this study before Al 
coating [29], [SOj. In the top panel of Fig. [4j we exam- 
ine (T) c>c Us (T) /m* as a function of temperature in 
four different samples; with x = 0.038 corresponding to 
the two far underdoped samples having different mea- 
sured values of Aab(O) in the region of coexisting phases, 
X = 0.074 being close to optimal doping and x — 0.10 be- 
ing an overdoped concentration, all of which were used 
to determine the values of Aaf,(0) shown in the top panel 
of Fig. [31 It should be noted that the orange and black 
curves in Fig. [4lfor x = 0.038 were made using the same 
AXab{T) data, but different values of Aab(O) because the 
temperature dependence of only one of the two samples 
shown in the top panel of Fig. [3l was measured before 
aluminum coating. As can be seen in the top panel of 
Fig. m the values of X~^{T 0) are quite large for 
X > 0.047 (red diamonds and green triangles) relative 
to those with x < 0.047 because the measured values of 
Aab(O) are much smaller than 1 /xm, i.e. 0.182 nm and 
0.270 nm. However, for x < 0.047 (black circles and 
orange squares) the values of A~^^(0) vary much less be- 
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FIG. 3: (Color online) Top panel: The zero temperature Lon- 
don penetration depth, Xab{0), as a function of the Co con- 
centration, X. The three dashed blue lines are theoretical 
curves obtained using Eq.[2]for three different values of XabiO) 
in the pure superconducting state. The solid gray line is a 
guide to the eye (in the form of A+B/x"). Also shown are 
values of Xab{0) obtained by other experiments for compari- 
son explained in the text. Bo ttom panel : Phase diagram for 
Ba(Fei_,Co,)2As2 [l,i,[l,Wr 



cause Aab(O) becomes closer to 1 /im, i.e. 0.673 nm and 
0.921 nm. 

Using the same penetration depth data that was used 
in the top panel of Fig. [31 we construct the normalized su- 
perfluid density (phase stiffness), psiT) = (A(0)/A(^))^ 
which is commonly used to analyze penetration depth 
data and a quantity which is fairly easy to calculate for 
an arbitrary gap structure. The bottom panel in Fig. [31 
shows ps (T) for the same samples shown in the top panel. 
The black and orange curves were constructed using the 
same AXab{T), but different Aab(O) values for the heav- 
ily underdoped sample, a;=0.038, and the red and green 
curves are the data for optimally doped and overdoped 
compositions, respectively. Also shown for comparison 
are the Ps{T) curves for a single band s-wave supercon- 
ductor (dotted blue line) and a d-wave superconductor 
(dotted gray line), both in the clean limit. From Fig. [H 
Ps{T 0) and ps{T Tc) behave quite differently for 
the members of the Ba(Fei_a;Co3;)2As2 series compared 
to the standard, single gap s-wave and d-wave clean limit 
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cases. Impurity scattering would turn the d-wave curve 
quadratic at low temperatures, while leaving s-wave al- 
most intact. 

The data for all doping levels show an overall similar 
trend of the evolution of ps (T) across the phase diagram. 
A special feature of this behavior is the negative curva- 
ture just below Tc- This behaivor suggests that below 
Tc the superconducting gap develops slower than it does 
in the case of a single gap, which is a feature consistent 
with the behavior of Ps(T) in a two-gap superconduc- 
tor [3^ . Furthermore, the normalized Ps(T) curve for 
the optimally doped sample over the entire temperature 
range stays above the curves for both heavily underdoped 
and overdoped samples. This distinction between the 
different Co-doping compositions suggests that the gap 
anisotropy, which is generally considered as being either 
the actual angular variation in fc— space and/or the de- 
velopment of an imbalance between the gaps on differ- 
ent sheets of the Fermi surface, notably increases in the 
overdoped and underdoped compositions. This is con- 
sistent with the measurements of the specific heat jump 
[33! and the residual term 34], as well as with measure- 



ments of thermal conductivity |35l . |36|. Thermal con- 
ductivity measurements with heat flow along the c-axis 
actually suggest that nodal regions develop in the super- 
conducting gap in heavily under- and over-doped compo- 
sitions. Indeed, measurements of Ac in a closely related 
Ba(Fei_a;Ni:c)2As2, also suggest the existence of nodes in 



•-3.8% Co (X„=0.921 urn) 
■ - 3.8% Co (X„=0.673 urn) 
• - 7.4% Co (^„=0.270 urn) 
•-10% Co (X„=0.182 urn) 



the superconducting gap [37 1 




CONCLUSION 

In conclusion, the zero temperature value of the in- 
plane London penetration depth, Aab(O), has been mea- 
sured for the Ba(Fei_j:Coj:)2As2 series across the super- 
conducting "dome" of the phase diagram using an Al 
coating technique along with TDR measurements. There 
is a clear increase in Aab(O) below x ~ 0.047, which is 
consistent with a reduction in the superfluid density due 
to the competition between itinerant antiferromagnetism 
and superconductivity for the same electrons. The mea- 
sured values of \ab (0) were also used to construct the nor- 
malized superfluid density (phase stiffness), Ps{T), and 
study its evolution with doping. The negative curvature 
of Ps{T) just below Tc for samples across the supercon- 
ducting dome of the phase diagram implies two-gap su- 
perconductivity. A notable suppression of ps for heavily 
underdoped and slightly overdoped samples with respect 
to samples with optimal doping suggests a developing 
anisotropy of the superconducting gap toward the edges 
of the superconducting dome, consistent with the behav- 
ior found in specific heat and thermal conductivity stud- 
ies. 
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FIG. 4: (Color online) Top panel: X~^{T) for samples of dif- 
ferent Co concentrations, x, from the Ba(Fei-2,Coa;)2As2 se- 
ries constructed from previously measured A\ab{T) curves 
and Aa!,(0) from this study. Bottom panel: Normalized su- 
perfluid density, Ps{T), for the same samples shown in the 
top panel along with the standard s-wave and d-wave cases 
for low impurity scattering. 
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